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ABSTRACT 

This paper reviews flat plate skin friction data 
from early correlations of drag on plates in water 
to measurements in the cryogenic environment of 
The NASA Langley National Transonic Facility 
(NTF) in late 1996. The fiat plate (zero pressure 
gradient with negligible surface curvature) 
incompressible skin friction at high Reynolds 
numbers is emphasized in this paper, due to its 
importance in assessing the accuracy of 
measurements, and as being important to the 
aerodynamics of large scale vehicles. A 
correlation of zero pressure gradient skin friction 
data minimizing extraneous effects between tests 
is often used as the first step in the calculation of 
skin friction in complex flows. Early data compiled 
by Schoenherr for a range of momentum thickness 
Reynolds numbers, R 01 from 860 to 370,000 
contained large scatter, but has proved 
surprisingly accurate in its correlated form. 
Subsequent measurements in wind tunnels under 
more carefully controlled conditions have provided 
inputs to this database, usually to a maximum Rg 
of about 40,000. Data on a large axisymmetric 
model in the NASA Langley National Transonic 
Facility extends the upper limit in incompressible 


R e to 619,800 using the van Driest transformation. 
Previous data, test techniques, and error sources 
are discussed, and the NTF data will be discussed 
in detail. The NTF Preston tube and Clauser 
inferred data accuracy is estimated to be within 
- 2 percent of a power-law curve fit, and falls 
above the Spalding theory by 1 percent at R 0 of 
about 600,000. 


1. INTRODUCTION 

The design of transport aircraft requires that 
accurate estimates of skin friction be made at 
length Reynolds numbers around 10 9 and Mach 
numbers of approximately 0.8, corresponding to 
cruise conditions 1 . This estimate is often made by 
first calculating the flat plate incompressible skin 
friction and then correcting for various effects such 
as pressure gradient, three-dimensionality of the 
flow, compressibility, etc. Several baseline 
theories/correlations are available, most of which 
differ in the skin friction level predicted at high 
Reynolds numbers where there has been a dearth 
of data. 

As will be shown, the most commonly used 
correlations of skin friction do not agree to the 
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desired accuracy at high Reynolds numbers, and it 
is not clear which method is the most accurate. 
Over the years, compilations and critical reviews 
of the available data have been made, but only 
recently has skin friction data at very high 
Reynolds numbers become available. The most 
recent data is usually tunnel wall data, with the 
exception of the present data which was taken on 
a large axisymmetric model. 

An experiment to measure flat plate skin friction 
at flight Reynolds numbers was conducted in the 
Langley National Transonic Facility, NTF. The 
purpose of the test was to provide skin friction 
data at very high Reynolds numbers to compare 
with existing theories. The difficulties in 
measuring skin friction free of extraneous effects 
and to the desired accuracy was recognized early. 
In addition to the usual problems associated with 
testing at high Reynolds numbers, the cryogenic 
environment of the tunnel affected both the model 
and instrumentation and complicated the test. A 
description of the test, the results obtained, and 
comparison with other available skin friction data 
will be discussed. 

A factor to be considered in testing at high unit 
Reynolds numbers is the difference in scales 
between the wall and the outer region. At the wall, 
the turbulent scales are so small that wall 
roughness may be an important factor — /.e., the 
wall must be very carefully machined to avoid 
roughness-induced effects. What for most 
engineering purposes would be considered 
adequate wall smoothness may not be the case at 
high Reynolds numbers 2 . Another of the many 
problems is the difficulty of manufacturing 
boundary layer rakes with tubes spaced closely 
enough to capture the law-of-the-wall region. 


2. NOMENCLATURE 
c f friction coefficient 

c p pressure coefficient 

C constant in equation 7 

d outside diameter of Preston tube 

D model diameter, 12.75 inches 

E constant in equation 2 

H shape factor, §*/0 

k van Driest constant, equation 7 

M Mach number 

r radial coordinate 

R unit Reynolds number, pu/p 


s distance along surface of model 
T temperature 

u velocity, fps 

u T shear velocity, (x w /p) 1/2 

u + u/u T 

x axial distance along model 

y coordinate normal to surface 

y* y uA> 

a angle of attack of model 

P Clauser pressure gradient parameter, 

(5/x ) (dp/dx) 

5 boundary layer thickness 

S u thickness at u /Ue =0.99 

5 displacement thickness 

5 + value of y + at edge of boundary layer 

0 momentum thickness, inches 

p dynamic viscosity 

V kinematic viscosity 

p density 

1 shear stress 

<(> roll angle on model, see Figure 4 

subscripts 

aw adiabatic wall condition 

e at edge of boundary layer 

i incompressible 

t stagnation condition 

w at wall 

x based on x coordinate 

0 based on 0 

freestream condition 

superscript 

reduced to incompressible form by Van 
Driest method (equation 6) for velocities and 
integral quantities, or by Sommer and Short T 
method (equation 5) for temperature-dependent 
quantities 

3. NTF EXPERIMENT 

The NTF test was designed to accurately 
measure adiabatic flat plate skin friction values at 
Reynolds numbers as high as possible and Mach 
numbers corresponding to flight conditions in order 
to extend the existing skin friction-Reynolds 
number database. The highest Reynolds numbers 
obtainable in Langley Research Center tunnels 
are produced in the cryogenic transonic tunnels, 
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NTF and 0.3 M Transonic Tunnel. For this test, 
unit Reynolds numbers as high as 94 x 10 6 /ft were 
run, and length Reynolds numbers of 940 x 10 6 at 
the downstream measurement station on the 
model were measured. The difficult task of 
measuring skin friction in the cryogenic tunnels is 
made even more difficult by model and 
instrumentation problems related to the extreme 
low temperature of the flow. 

The idea of a test for obtaining high Reynolds 
number Cf data in NTF is not new, and in fact, a 
program to do so was outlined by Saric 3 . The 
ideas presented in this report were used in 
planning the NTF test, the major difference being 
the model on which Cf was to be measured. It was 
determined in the early planning stage that a two- 
dimensional flat plate posed too many problems in 
mounting and maintaining the accuracy of the 
surface in the high dynamic pressure environment. 
For this reason an axisymmetric model was 
designed for which transverse curvature effects 
were small, based on the thickness of the 
predicted boundary layer. In addition, it was 
necessary to account for compressibility effects at 
Mach numbers as high as 0.8 in order to 
transform the data to an equivalent incompressible 
state. 

It was decided that no new skin friction 
measuring techniques would be developed for this 
test, since development was considered too time- 
consuming for an experiment in NTF. It was also 
required to measure the skin friction as accurately 
as possible, a notoriously difficult measurement to 
make, even at ambient temperature. Three 
different standard techniques were used with the 
rationale that they would either provide 
consistency among themselves or point out 
inadequacies in the experiment. A skin friction 
balance was used since it is the only way to 
measure skin friction directly. Also used were 
Preston tubes and boundary layer surveys, from 
which skin friction was inferred by a modified 
Clauser method. There has been a renewed 
interest in the validity of the law of the wall lately, 
as evidenced in the report of George and Castillo 4 . 
The law of the wall is used extensively here to 


present data and to infer Cf from velocity profile 
measurements. The surveys were also necessary 
to determine the boundary layer integral quantities 
such as the momentum and displacement 
thicknesses. 

A photograph of the model in NTF is shown in 
Figure 1. 

4. FLAT PLATE SKIN FRICTION DATA AND 
THEORY 

4.1 Theories and Correlations 

Schoenherr s correlation of his own data and 
the data of others dating to 1872 5 was published 
almost 70 years ago, but is still used in its 
correlated form, except at low Reynolds numbers. 
Skin friction was obtained from the drag of plates 
of various sizes towed in water and, not 
surprisingly, the data exhibits large scatter when 
plotted as a function of Reynolds number. When 
correlated with an equation of von Karman, the 
results agree with most carefully-controlled 
experiments. Twenty years later, in 1953, the 
paper by Landweber was published 6 , describing 
the characteristics of turbulent boundary layers 
and deriving a shear law that agreed closely with 
the method of Schoenherr. Nine years after this, 
the method of Spalding 7 was published. This 
method was derived from his sublayer-buffer-log 
profile, but did not account for the outer portion of 
the boundary layer. If the usual law of the wall 
constants are used in Spalding s Cf equation the 
answer is incorrect, since the wake is not 
accounted for correctly. The constants were 
adjusted for use as the incompressible theory in 
the Spalding-Chi compressible skin friction 
method 6 to produce a more reasonable Cf level. 
By adjusting the constants of course, the level 
can be changed, but, more important, the slope of 
Cf vs R fi does not match that of the Karman- 
Schoenherr correlation over the complete 
Reynolds number range. Other methods which 
have been used are Ludwieg-Tillmann 9 , and 
recently, Fernholz and Finley 10 , which is based on 
both inner and outer similarity and is compared 
with data to R 0 of 200-300 x 10 6 . 


The Karman-Schoenherr equations are: 

0.242 / = Io9io(cf^x) ~ ^9io(^^e) 

c f =0.558c f /| 0.558 + 2^/c^j 


( 1 ) 
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where Cf is the average skin friction coefficient. 

The local skin friction coefficient, Cf, was obtained by differentiating the first of the two equations. 

The Spalding equations are 


R = 
x 12 
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( 2 ) 


The Ludwieg-Tillmann equation is 

c, = 0.246 x 10-°- 678H x R-° 268 


(3) 


The Fernholz and Finley equations are 


A- Uln 
2 Ik 


where In 


AUt ' + C-N + I — — M llnf ^ 


^ Au A 


= ln(R x ) = 0.3 + ln(R e ) 


A = f y ^2 dy 

Jo u 


and In — =-2.7 for R 0 > 2000 


= -0.404 ln(R e ) + 0.37 for 4252 < R e < 2000 
k = 0.4, C = 5.1, M = 4.70, and N = 6.74 


(4) 


The methods of Karman-Schoenherr, 
Spalding, Ludwieg-Tillmann, and Fernholz and 
Finley are shown in Figure 2. From the figure, it is 
obvious that the theories do not agree with each 
other over the complete range of Reynolds 
numbers. The methods of Karman-Schoenherr 
and Spalding show opposite trends at low and 


high Reynolds numbers, the crossover point 
occurring at R 0 between 6000 and 7000. The 
Karman-Schoenherr correlation includes Kempf s 
data to R e of about 370,000 (R x of 450 x 10 6 .) 
The Spalding method relies on the constants of 
the law-of-the-wall, and thus might be in question 
at high Reynolds numbers unless these constants 
are shown to be the same as for lower Reynolds 
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numbers. Ludwieg-Tillmann deviates from 
Karman Schoenherr below R 0 about 3000 and 
above R 0 about 20000. Fernholz and Finley is 
lower than the other theories. 

4.2 Data 

Over the past 30 years excellent reviews of 
available incompressible and compressible 
boundary layer data have been compiled. These 
include the compilation for incompressible flows by 
Coles and Hirst from the Stanford Conference in 
1968 11 , the compilations of compressible flow data 
by Fernholz and Finley in 1977 12 , and 1981 13 , and 
the examination of Incompressible data to 1996 by 
Fernholz and Finley 10 . Since that time, other data 
has been published, most high Reynolds number 
data having been taken on the wall of wind 
tunnels. These data appear to be different from 
flat plate boundary layer data 14 , i.e., the boundary 
layer has developed in regions of wall curvature 
and strong adverse and favorable pressure 
gradients, and is usually not characteristic of 
equilibrium flat plate boundary layers until a 
sufficiently long run has been made. Relaxation 
has apparently not been completed before 
entering the test section. 

Reference 15 from 1996, listed data in 
addition to that of reference 10, but did not present 
skin friction measurements at high Reynolds 
numbers. In 1984, Gaudet 16 , published data at 
Mach 0.8 on the sidewall of the RAE 8 ft by 8 ft 
wind tunnel and reexamined the data of reference 
17. Other data was published in 1994 by 
Motallebi 18 on the wall of a wind tunnel for R e from 
26,000 to 106,000 in compressible flow; however, 
skin friction was not measured in ail cases and the 
data are not used for comparison here. The 
tabulated data of Fernholz, Krause, Nockemann, 
and Schober 19 is compared in Figure 3 with high 
Reynolds number data of Gaudet, low Reynolds 
number data of Coles 20 and Purtell 21 , and the 
methods of Karman-Schoenherr (K-S) and 
Spalding. For this figure the compressible profiles 
of Gaudet, were re-reduced and transformed by 
the van Driest transformation 22 , which was also 
used for the NTF data transformation. Noted in 
the figure is the R e range of the data taken in NTF 
on an axisymmetric model, and reported in a 
following section. 

The methods of Karman-Schoenherr and 
Spalding will be used to compare with the present 
data, since they agree reasonably well with the 
best experimental data available at this time over 
a large range of FV 


5. NTF EXPERIMENT -RESULTS AND 
DISCUSSION 

5.1 Tunnel. Model. Instrumentation, and Test 
Conditions 

5.1.1 Tunne l The range of operating 
conditions for the National Transonic Facility 
(NTF) at Langley Research Center 23 is as follows: 
Mach numbers from 0.2 to 1.2, total pressure from 
one to 9 atmospheres, total temperature from 
— 320 to 150deg. F, and Re/ft from 3.7 to 146 x 
10 6 at Mach 1. For cryogenic operation, liquid 
nitrogen is injected into the flow downstream of the 
test section, and vaporized to maintain low tunnel 
temperatures. Intermediate temperatures can be 
attained by regulating the amount of liquid nitrogen 
injected. The tunnel can also be run using air if no 
cooling is required. For the present tests, both 
ambient temperatures and cryogenic temperatures 
were used to obtain the desired Reynolds number 
range. 

The test section is 25 feet long, and 8.2 ft 
square 24 , making it possible to mount very long 
models in the tunnel to produce large length 
Reynolds numbers. The model described here 
was 17.28 feet long. The test section is slotted to 
minimize blockage effects, which were 3.3% 
based on an inviscid geometric area ratio. 
Measured model pressures demonstrated that the 
flow over the model was uniform, validating that 
the slotted design worked as expected. 

Standard tunnel data reduction and tunnel 
instrumentation 25 were used; however, local 
values of flow parameters on the model were 
recalculated from measured model conditions 
using the Fortran routines of NTF. The Beattie- 
Bridgemann equation of state was used to 
calculate the properties of both nitrogen and air 26 . 

5.1.2 Model The design of the model was a 
compromise among several factors. It was 
desirable to have a long model in order to produce 
large length Reynolds numbers. The diameter 
was required to be large enough so that the model 
would not deflect significantly under its own weight 
and to allow room for the balance and other 
instrumentation to be mounted internally. In 
addition, it needed to be large enough that 
transverse curvature effects would not significantly 
affect skin friction. It could not be too large, or 
model loads and blockage effects would be 
excessive. 

A sketch of the axisymmetric 347 stainless 
steel model is shown in Figure 4. It was an 
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axisymmetric cylinder 12.75 inches in diameter 
having a nose described by a superellipse, and 7 
cylindrical sections downstream of the nose. Ports 
were installed in the model at Stations 1 and 2, at 
x=73.95 and 121.95 inches, respectively. Skin 
friction was measured by a balance and by 
Preston tubes at Station 2, and by rakes at 
Stations 1 and 2. The whole model, including the 
nose, was polished to a surface finish of 4 p in. 
No transition trips were used on the nose. The 
model was sting-mounted in order to minimize 
boundary layer interference effects, which would 
be caused by model struts. The NTF sting mount 
had roll and angle of attack capability that could be 
used to set the model at nominally zero angle of 
attack. Initially, Preston tubes were mounted 
circumferentially around the model to measure the 
nonuniformity of the flow and adjust the angle of 
attack. 

The design of the nose and flow over the 
model were predicted using CFD methods. It was 
determined from these calculations and boundary 
layer calculations 27 that the pressure gradient 
would not be a factor in measuring flat plate skin 
friction. Similarly, the ratio of the boundary layer 
to the model radius was estimated 27 to be 0.25 at 
the second measurement station, and transverse 
curvature effects on Cf were estimated to be less 
than 1.5 %, at R e of about 600,000. 

5.1.3 Instrumentatio n Three methods of 
determining skin friction on the model were used: 
a skin friction balance, Preston tubes, and velocity 
profiles from which skin friction was inferred by the 
Clauser method. Photographs of the three 
instruments used on the model are shown in 
Figure 5. Preston tubes and the skin friction 
balance were tested at Station 2 at both hot and 
cold flow conditions. The boundary layer rake was 
used at Station 1 at the hot flow condition and at 
Station 2 at hot and cold flow conditions. 

The problems associated with the direct 
measurement of skin friction by balance are well 
known . Additional problems resulting from 
operating a mechanical device in a cryogenic 
environment required that a risk reduction 
experiment be done in a smaller tunnel — the 
Langley 0.3 Meter Transonic Cryogenic Tunnel. 
The problems encountered in the risk reduction 
experiment were larger than anticipated and 
greatly influenced the design of the final balance. 

Preston tubes were sized from boundary layer 
calculations. They were designed so that the tube 
was in the logarithmic region of the boundary layer 


at all test conditions. It has been shown 29 that the 
tube can extend well into the wake portion of the 
boundary layer with no adverse effects; however, 
the tubes were 0.058 inches in diameter, and did 
not extend into the outer portion. The nose was 
unchamfered and polished, and the ratio of 
internal to external diameter was 0.6. Care was 
taken in mounting the probes so that the tubes 
would be firmly attached to the surface and to 
ensure that no mounting apparatus would protrude 
into the flowfield (see Figure 5 b). 

The only set of Preston tube data know to the 
authors covering the range of Reynolds numbers 
for which measurements were to be made was 
that from Laval University 30 . Rather than rely on a 
single set of data, a high Reynolds number 
calibration was made at Princeton University in the 
Superpipe Facility 2 . The data from the Princeton 
test was cast in the parameters of four different 
calibration methods 31 ' 4 to see if any one method 
gave superior results in determining skin friction 
from known inputs. It was found that all methods 
gave the same results. Patel s high Reynolds 
number calibration was found to predict the 
Princeton data at high Reynolds numbers, and 
was used to reduce the NTF data. 

Inference of skin friction from velocity surveys 
is most accurate if tubes on the boundary layer 
rake are positioned within the logarithmic portion 
of the turbulent boundary layer. In addition, the 
edge of the boundary layer should be defined so 
that the integral quantities 6 and 6 are accurately 
calculated from the profiles. The boundary layer 
calculations discussed previously provided the 
guidelines for tube placement. The tips of the pitot 
tubes were not chamfered, and were polished, as 
was done for the Preston tubes 

Model pressures were measured by an 
electronically scanned system, and model 
temperatures by copper-constantan (Type T) 
thermocouples accurate to 1 degree K. One port 
on each ESP module was dedicated to measuring 
a known reference pressure that was also 
measured by a secondary standard. Whenever 
the reference pressure difference deviated by 
- 0.19% of full scale the modules were re-zeroed 
online. 

5.1.4 Test Conditions The model was tested at 
Mach numbers from 0.2 to 0.85 and unit Reynolds 
numbers from 6 x 10 6 to 94 x 10 6 per foot. A 
matrix of the test conditions is shown in Figure 6. 
The highest unit Reynolds number was attained at 
Mach 0.6. The highest unit Reynolds number at 
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which data could be obtained at Mach 0.85 was 65 
x 10 6 , a consequence of the load constraints for 
the model. Noted in the figure are the conditions 
at which data was obtained on Preston tubes, the 
balance, and the boundary layer rake. At most 
tunnel conditions, all three devices were tested, 
but at Mach 0.7 only Preston tubes were run. 

Due to the finite blockage of the model in the 
tunnel, it was necessary to calculate the local 
Mach number, which could be slightly different 
from the freestream Mach number. This was done 
by averaging pressures at three orifices just ahead 
of Station 1 and using this pressure as the local 
static pressure to calculate the nominal local Mach 
number. To define the model pressure in the 
vicinity of the measurement stations, 11 orifices 
around Station 1 and 11 around Station 2 were 
averaged. These averaged pressures were used 
as the local static pressure in calculating c p at their 
respective stations. Similarity, three temperatures 
near the two primary data positions were averaged 
in calculating Twl and Tw2. 

Representative pressures and temperatures 
over the model are shown in Figure 7. It is 
evident that the pressure gradient is sufficiently 
close to zero and that its effect on the local shear 
stress can be considered negligible. This is more 
accurately quantified in later discussion. It is also 
evident that there is little change in the pressure 
gradient with Mach number. While not as uniform 
as the pressure distribution, the model 
temperature variation is small, being more uniform 
at ambient temperature, as would be expected. 

5.2 Factors affecting NTF Skin Friction 
Measurements 

There are at least seven factors associated 
with the conditions of the NTF test that can have 
an effect on measured surface shear. These 
factors are compressibility, pressure gradient, heat 
transfer, surface curvature, surface roughness, 
lateral divergence (three-dimensionality) of the 
flow, and the freestream disturbance level and 
scale. All effects, except the freestream 
disturbance level, were examined before the NTF 
test was begun to insure that they could be 
controlled/minimized to the extent that a flat plate 
skin friction could be measured. 

5.2.1 Compressibility Effects T ests were 

conducted at Mach numbers from 0.2 to 0.85. 
Above approximately Mach 0.3, compressibility 
effects must be considered in the flow. One 
approach, and the one used by Allen in Preston 
tube reductions 29 31 , is the use of the T method. In 


this method temperature-dependent quantities 
such as p and p are replaced in parameters and 
equations by their values at a temperature 
intermediate between the wall and freestream 
values -- the T-prime condition. Allen used the 
Sommer and Short method, which is the Rubesin 
and Johnson equation with constants determined 
for Mach 3.82 in air 35 . The Sommer and Short T- 
prime method was used in this report for all 
Preston tube reductions. The equation defining 
the T condition is 


T. = 1 . + 0.035 M? + 0.45 
T 

' p 



(5) 


In order to present velocity profiles and 
integral quantities in an equivalent incompressible 
form, velocity profiles were transformed using the 
Van Driest method 21 , which was shown to 
transform velocity profiles in the outer region of 
adiabatic wall boundary layers at Mach 11 in 
helium 36 . Velocities are transformed in proportion 
to the density variation from the wall to the edge of 
the boundary layer by this method; i.e., the vertical 
coordinate, y, is not transformed. 

Called generalized velocities by Maise and 
McDonald 37 , the van Driest transformation is 
defined as 



With the transformed velocity profile and the 
density evaluated at the wall, new values of 6,8* 
and u e result. These transformed properties are 
used to calculate R 0 and c f . The transformation 
was coupled with the Clauser method to infer the 
shear stress from velocity profiles by an interactive 
graphical method. The incompressible profile was 
plotted along with the incompressible law of the 
wall, and the data was fit to the law of the wall by 
iteration. It is estimated that the data could be fit 
to the law of the wall with a resolution of about 
0.5 % in Cf by this method. The accuracy of the 
results are a function of the constants in the law of 
the wall, defined as 


where k=0.41 and C=5.5 

The profile data were originally reduced with 
C=5.0, however it was found that C=5.5 gave a 
better fit to the Preston tube data and skin friction 
theory. 


V u , 


+ C 


(7) 
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5.2.2 Pressure Gradient-Relaxation Effects 
Preliminary tests were run in the Langley 0.3 
Meter Transonic Cryogenic Tunnel, where the 
effects of both favorable and adverse pressure 
gradients on measured skin friction were 
examined. This was accomplished by moving the 
upper wall of the test section to produce pressure 
gradients of different magnitude. Figure 8 (a) 
shows the effects on Preston tube measured Cf in 
terms of (3, the pressure gradient parameter used 
by Clauser to generate self-similar equilibrium 
profiles 38 . Data are shown at Mach numbers of 
0.6 and 0.8 at two different unit Reynolds 
numbers, and pj is a transformed value obtained 
by using the van Driest transformed value of 5 . 
Compressibility and Reynolds number effects are 
not completely removed using p> however, the 
correlation should be adequate for estimating local 
pressure gradient effects. For these test 
conditions 

— = 1-0.7 ft (8) 

fi-o 

As shown in Figure 7(a), the pressure gradient 
on the nose is strongly favorable switching to 
strongly adverse about x=13 inches, and relaxing 
to very small adverse pressure gradient about 
x=30 inches. The pressure gradient is zero at 
x=50 inches. The question arises as to whether 
the outer part of the boundary layer had relaxed to 
the equilibrium zero pressure gradient state by the 
first measurement station. Station 1 was located 
at x=73.95 inches, 44 inches from the beginning of 
small adverse gradient. Measured values of 8 and 
0 were 1.3-15 and .065-. 075 inches here, or 31 8 
and 630 0 from the beginning of small adverse 
pressure gradient and 17 8 and 340 0 from zero 
pressure gradient. The distance required for 
relaxation for flows with adverse pressure gradient 
followed by zero pressure gradient is about 5-10 
5 39 

5.2.3 Heat Transfer Effects The effect of 
nonadiabatic wall temperature effects on skin 
friction were also measured in the 0.3 Meter 
Tunnel. The total temperature of this tunnel can 
be rapidly varied from above ambient temperature 
down to 100 K, making it possible to measure Cf 
over a range of T w /T aw The results are shown in 
Figure 8 (b) for Mach 0.6 and Mach 0.8. 
Repeatability is shown by the agreement between 
two different runs. For most runs in this tunnel, 
the wall was slightly nonadiabatic due to heat 
transfer through the tunnel wall to a plenum 


chamber behind the wall, and then to outside the 
tunnel. 

5.2.4 Transverse Curvature Effects The 
transverse curvature effect is a geometric effect 
characterized by the ratio of 8 to model radius. 
Boundary layer calculations were made using 2- 
foot long elliptical noses, model pressure 
distributions from Newtonian impact theory, 
estimated transition locations, and model radii of 3 
and 6 inches to estimate the magnitude of this 
effect. The measured ratios of dr were found to 
be 0.25 at the second data acquisition station, 
resulting in an increase in Cf less than or equal to 
15%, based on the calculations. 

5.2.5 Surface Rouahness Effects In order to 
insure that the model was hydraulically smooth, 
boundary layer calculations were used to 
determine the y + = 5 height normal to the surface. 
Figure 9 shows y in inches plotted against y + for 
the most severe case, that at highest unit 
Reynolds number. Based on this figure, the rms 
wall roughness should be smaller than 30 pin. 
When constructed, the complete model and skin 
friction balance element had a measured 
roughness of 4 p inches and a value of y + of about 
one. Based on the y + =5 criterion, surface 
roughness should not have significantly affected 
the skin friction measurements of this test. 

5.2.6 Lateral Divergence Effects Three 
dimensionality of the flow is evident even in flat 
plate experiments where great care has been 
taken to insure the two-dimensionality of the outer 
flow. From velocity profiles measured across the 
flow, plots of 0 are always somewhat ragged, in 
keeping with skin friction measurements at the 
same locations and variations in the wake 
component. This may be due to the raggedness 
in the transition process, and/or the presence of 
near-wall structures in the boundary layer, 
inducing a standing pattern in the flow. When 
there is finite three dimensionality in the flow, 
effects are added to the usual two dimensional 
raggedness. In order to minimize the three 
dimensionality of the flow on a longitudinal 
cylinder, the expected droop due to the model 
weight was calculated and found to be 
approximately 0.25 degree. Since the model 
could be adjusted within 0.1 degree, it was felt 
that the effect could be minimized by careful 
alignment of the model. Initial runs were devoted 
to model alignment. In addition, rake, balance, 
and Preston tube measurements were always 
made in the same azimuthal plane, which will be 
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shown to be an important issue in a later section 
of this report. 

5.2.7 Freestream Disturbance Levels The 
dependence of turbulent Cf and other boundary 
layer parameters on rms freestream turbulence 
intensity and length scale has been shown to be 
nonlinear. 40 For NTF the freestream turbulence 
level is not known, so this parameter remains a 
potentially important parameter in the assessment 
of skin friction measurements made in this tunnel. 
The main sources of noise in NTF have been 
identified 41 by measuring surface pressure 
fluctuations at several places within the tunnel. It 
is noted that, in comparison with other transonic 
wind tunnels, the NTF has low levels of test 
section fluctuating static pressure. 

5.3 Profile Measurements — Skin Friction from 
Clauser Method 

5.3.1 Transformed Velocity Profiles The pitot 
rake was mounted at both Stations 1 and 2. Skin 
friction coefficients were inferred from velocity 
profiles by first transforming the data to 
incompressible form by applying the van Driest 
transformation and then using the Clauser method 
of fitting the velocity profile to the law of the wall. 
Pitot pressures were reduced to Mach numbers 
assuming isentropic, ideal gas flow, constant 
gamma, and constant p s throughout the boundary 
layer 42 . Since T t was not measured, the local total 
temperature in the boundary layer was estimated 
by the Rotta relation 43 . By this means, both 
temperature and velocity profiles can be obtained 
from the Mach number profile along with 
measurement of T t and T w . Despite efforts to 
locate rake tubes very near the model wall, the 
first point in the velocity profile was well within the 
logarithmic layer of the boundary layer, as shown 
in the velocity profile data of Figure 10 (a) plotted 
in law-of-the-wall coordinates. In order to 
accurately calculate boundary layer integral 
quantities such as 5* and 0, the theoretical profile 
of East, also used by Gaudet 16 , was used to fill in 
the profile from the wall to the first data point. The 
profile of Spalding 7 could also have been used. 
For the law of the wall k=0.41 and C=5.5 were 
used. The strength of the wake component was 
2.2 for the present data, in agreement with other 
high Reynolds number data as shown in 
Reference 44. Transformed profiles are plotted in 
traditional defect coordinates in Figure 10 (b) and 
in the coordinates of Reference 4 in Figure 10 (c). 


Transformed shape factors are shown in Figure 
11 along with the low Reynolds number data of 
Coles and Purtell, the higher Reynolds number 
data of Gaudet, and the results from a zero 
pressure gradient boundary layer calculation. For 
the calculation of shape factors for data, the value 
of C used did not affect the values of transformed 
shape factors significantly. 

5.3.2 Inferred Skin Friction Figure 12 shows 
values of Cf inferred from velocity profiles by the 
Clauser method plotted as a function of 
transformed R e . The value of skin friction 
obtained is dependent on the constants used in 
equation 7. The value of k is usually taken to be 
0.41 in the literature and was used for reducing 
the present data. The value of C chosen by 
various authors to describe their data usually 
varies between 5.0 and 5.5. This variation is not 
surprising since, as pointed out by Coles, the outer 
part of the boundary layer is extremely sensitive to 
extraneous effects in the flow. A value of C=5.5 
was used for the final reduction of the profile data 
since this value gave better agreement with 
Preston tube data, the theory of Spalding, and the 
results from finite-difference calculations than did 
5.0. To decrease the Cf level to that of Karman- 
Schoenherr, an even higher value of C would have 
been necessary. Reduction of the data using 
C=5.0 and 5.5 is shown in Figure 12. A value of 
5.0 gave values of Cf too high at Reynolds 
numbers between 30,000 and 100,000 where 
other data exists, and for this reason a value of 5.5 
was used for the final reduction of the present 
data. 

5.3.2 Circumferential Measurements To check 
for non-uniform flow effects, rake data were taken 
at nominal a=0 at <(> increments of 15 degrees by 
rolling the model. Data from the pitot rake were 
reduced to boundary layer parameters, and are 
shown in Figure 13 (a). Measured skin friction 
coefficients and values of R e show a consistent 
behavior, /.e., as 9 decreases, Q increases; 
conversely, as 9 increases, Cf decreases. 

The Cf data are replotted as a function of R e in 
Figure13 (b). From this figure, it was concluded 
that when the skin friction is plotted as a function 
of R e , variations in circumferential flow are 
removed, and also that data from different 
measuring devices should be taken at the same 
circumferential location on a model. Two 
dimensional standing patterns on tunnel walls and 
models have been found, and the variation is not 
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always insignificant. For this test, skin friction data 
was always taken at 4>=0 and nominal ot=0. 

There is a significant gradient in the 
circumferential skin friction coefficient in the 
vicinity of <(>=0, where skin friction measurements 
were made. The Preston tube and rake 
measurements were essentially point 
measurements, however the skin friction balance, 
having an element 3 inches in diameter, 
subtended an angle of approximately 27 degrees. 
The skin friction was higher on the counter- 
clockwise side than at $=0, and lower on the 
clockwise side. 

To provide data for comparison with 
computational methods, additional runs were 
made with the model at angle of attack. Data are 
available from the authors. 

5.4 Preston Tube Data 

5.4.1 Model Alignment In order to align the 
model with the incoming flow, Preston tubes were 
mounted on plugs at six circumferential positions 
at the two skin friction measurement stations. 
Preston tube pressures, a measure of the local 
surface shear, have been shown to be more 
sensitive to flow angularity than surface 
pressures 18 . By reducing Preston tube data to 
skin friction values and plotting the results as a 
function of angle of attack, the model could be 
zeroed in relation to the flow, not the test section 
geometry. No yaw adjustment was possible with 
the tunnel sting-mounting arrangement; however, 
the model could be adjusted in angle of attack. 

Several runs were devoted initially to aligning 
the model, with results at Station 2 shown in 
Figure 14. It was determined that the nominal 
zero angle of attack was — 0.07 degrees, 
independent of tunnel Mach number. Later runs 
with the pitot rake installed and the model rolled in 
15 degree increments at a= — 0.07deg, confirmed 
these results by giving essentially the same 
pressure patterns on rake tubes 1 and 25 as the 
model was rolled. 

5.4.2 Skin Friction Measurements Preston 
tubes were mounted at Station 2 and data were 
taken in both ambient temperature air and cold 
nitrogen flow. The calibration of Patel was used, 
which was determined from the Princeton 
Superpipe calibration to be valid at the high 
Reynolds numbers of this test. Compressibility 
effects were accounted for by the T method of 
Sommer and Short 35 . 


Figure 15 (a) shows the resulting calculated 
compressible Cf at Mach numbers of 0.4, 0.6, 0.7 
and 0.85 as a function of tunnel Re/ft. The spread 
in Cf at the lowest unit Reynolds number between 
the Mach 0.4 and Mach 0.85 data is an indication 
of compressibility effects. In order to plot the data 
as a function of R 0) values of 0 at each test 
condition were calculated from the boundary layer 
survey data. The data in incompressible form are 
compared with the Karman-Schoenherr and 
Spalding theories in Figure 15 (b). 

5.5 Balance Measurements 

The skin friction balance was tested at Station 2 
at both hot and cold flow conditions. The resulting 
data are shown in Figure 16. The data are plotted 
as compressible Cf against tunnel unit Reynolds 
number in Figure 16 (a), and also as 
incompressible Cf against transformed R 0 in Figure 
16 (b) for comparison with theory. In 
incompressible form, the compressible Cf was 
transformed to incompressible by the T-prime 
method, and R 0 was taken from velocity profiles at 
the same location. In most cases the data appear 
to be significantly higher than incompressible 
theory, even at Reynolds numbers of 40,000 to 
50,000. The reasons for this could not be 
determined, however, in view of the severe 
problems encountered in running the balance at 
cryogenic temperature in the 0.3 Meter TCT, the 
same problems might be expected in NTF Also, 
the balance element footprint was much larger 
than that of the other devices, making it more 
susceptible to variations in the circumferential 
flow. 

5.6 Comparison of Three Methods 

Data from Preston tubes, velocity profiles, and 
the skin friction balance are shown in Figure 17. 
Data from the Preston tubes and velocity profiles 
are in good agreement. Balance data agrees with 
the other data at the highest Reynolds numbers, 
however there is large scatter in the balance data 
at R 0 of 40,000 to 200,000. 

The Preston tube and Clauser-inferred values of 
Cf agreed well with each other and were fit to a 
power law, having the following equation: 

c, = 0.0097 R e <_0144> ( 9 ) 

The data scatter about this line is within - 2%, 
as shown in Figure 18. Equation 9 is 1% above 
the Spalding value and 3% above the Karman- 
Schoenherr value at R e =600,000. It is equal to 
the Spalding value at R e =30,000. 
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This equation is for incompressible skin 
friction, being derived from the transformed NTF 
data in the range 32,680 < R e < 619,800. For 
application to compressible flow, a suitable T- 
prime equation could be used 35 . 

6. CONCLUSIONS 

Skin friction data at transformed R e of 619,800 
have been shown to agree with existing 
incompressible theories within 3% at the highest 
Reynolds number. At R e =30,000, the data are in 
agreement with the theory of Spalding and slightly 
above the Karman-Schoenherr correlation. The 
data obtained by two different methods were in 
good agreement; however, these data did not 
agree as well with direct measurements made with 
a skin friction balance. The scatter in the balance 
measurements, attributed to mechanical problems 
with the balance, is a strong factor in 
deemphasizing this data. 

The data is free from pressure gradient effects, 
but has transverse curvature effects which can 
increase Cf by as much as 1.5%. Surface 
roughness should not be an issue. A factor 
affecting the data may be the freestream 
turbulence of the NTF tunnel, both in rms 
turbulence level and in scale. There are 
indications that the turbulence level is not large, 
however its effect on the data cannot be 
quantified. 

The data were obtained from nominally 
incompressible flow at M=0.2 to compressible flow 
at M=0.85. The van Driest transformation and T- 
prime method appear to be adequate to remove 
the compressibility effects. 
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Figure 1 . Photograph of Model in NTF. 
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Figure 2. Flat Plate Skin Friction Theories 
and Correlations. 


Figure 3. Comparison of Flat Plate Data 
with Theory. 
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Figure 4. Sketch of Model. 




(a) Rake. 


Figure 5. Continued, (b) Preston Tube. 



(c) Skin Friction Balance. 


Figure 5. NTF Model Instrumentation. 
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Figure 6. Mach-Reynolds Number Test Conditions. 
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(a) Pressures. (b) Temperatures. 

Figure 7. Representative Pressures and Temperatures On Model. 
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Figure 8. Effects of Pressure Gradient and Wall Temperature on Measured Cf. 
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Figure 9. Calculated Turbulent Wall Scaling. 
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Figure 10. Transformed Velocity Profiles. 
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Figure 10. Transformed Velocity Profiles. 
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Figure 11. Transformed Shape Factors 


Figure 12. Effect of C on c* Obtained by 
Clauser Method. 
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Figure 13. Variation of Boundary Layer Properties from Rake with Roll Angle. 
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Figure 14. Determination of Nominal 
Zero Angle of Attack. 
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(a)Compressible Data. 
Figure 15. Preston Tube Data. 
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(b) Data in Transformed Parameters. 
Figure 15. Concluded. 


(a) Compressible Data. 
Figure 16. Balance Data. 
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Figure 18. fit to Preston Tube and Crofile Cf Data. 
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